Multi-channel pulse dynamics in a stabilized 
Ginzburg-Landau system 



H.E. Nistazakis 1 , D.J. Frantzeskakis 1 , J. Atai 2 , B.A. Malomed 3 , N. Efremidis 4 , and K. 

Hizanidis 4 

^Department of Physics, University of Athens, Panepistimiopolis, 15784 Athens, Greece 
School of Electrical and Information Engineering, The University of Sydney, Sydney, NSW 

2006, Australia 

^Department of Interdisciplinary Studies, Faculty of Engineering, Tel Aviv University, Tel 

Aviv 69978, Israel 

(^Department of Electrical and Computer Engineering, National Technical University of 

Athens, 15773 Athens, Greece 



Abstract 

We study the stability and interactions of chirped solitary pulses in a system of nonlinearly 
coupled cubic Ginzburg-Landau (CGL) equations with a group- velocity mismatch between them, 
where each CGL equation is stabilized by linearly coupling it to an additional linear dissipative 
equation. In the context of nonlinear fiber optics, the model describes transmission and collisions 
of pulses at different wavelengths in a dual-core fiber, in which the active core is furnished with 
bandwidth-limited gain, while the other, passive (lossy) one is necessary for stabilization of 
the solitary pulses. Complete and incomplete collisions of pulses in two channels in the cases 
of anomalous and normal dispersion in the active core are analyzed by means of perturbation 
theory and direct numerical simulations. It is demonstrated that the model may readily support 
fully stable pulses whose collisions are quasi-elastic, provided that the group- velocity difference 
between the two channels exceeds a critical value. In the case of quasi-elastic collisions, the 
temporal shift of pulses, predicted by the analytical approach, is in semi-quantitative agrement 
with direct numerical results in the case of anomalous dispersion (in the opposite case, the 
perturbation theory does not apply). We also consider a simultaneous collision between pulses 
in three channels, concluding that this collision remains quasi-elastic, and the pulses remain 
completely stable. Thus, the model may be a starting point for the design of a stabilized 
wavelength-division-multiplexed (WDM) transmission system. 
PACS: 42.81.Dp, 42.65.Tg, 42.81.Qb 
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1 Introduction 



It is commonly known that complex cubic Ginzburg-Landau (CGL) equations constitute a class 
of universal models for the description of pattern formation in various nonlinear dissipative media 
0| . Equations of the CGL type are also frequently used in nonlinear fiber optics, to describe 
the formation, stability, and interactions of solitary pulses (SPs). CGL equations with constant 
coefficients apply to a long nonlinear optic-fiber link if the pulses in it are broad enough, so that 
the corresponding dispersion length is essentially larger than the amplification spacing. In this 
case, the periodically placed amplifiers and filters (which limit the gain to a relatively narrow 
spectral band) may be considered in the uniformly-distributed approximation, neglecting their 
discreteness j2j. 

The single-component CGL equation with the cubic nonlinearity possesses a well-known ex- 
act solitary-pulse solution , which includes an internal chirp (phase curvature across the pulse) . 
However, a fundamental drawback of this solution is that it is unstable, as the zero solution to 
the cubic CGL equation, i.e., a background on top of which the pulse is built, is unstable due to 
the presence of linear gain in the equation. Development of physically realistic models in which 
solitary pulses are fully stable is a problem of an obvious interest in its own right, and it also has 
profound importance for fiber-optic communications (see the book and, for instance, a recent 
work as well as for the design of soliton-generating fiber-loop lasers ||. In the context 
of optical telecommunications, an issue of fundamental significance is the development of ade- 
quate models for multi-component systems, corresponding to a wavelength-division-multiplexed 
(WDM) multi-channel scheme implemented in the optical fiber. An objective is then to design a 
system supporting pulses in all the channels, which must be stable against small perturbations 
and mutual collisions (see, e.g., books |J, [f?J and a recent paper ||). 

A single-channel system, which may suppress the instability of the zero solution, simulta- 
neously allowing for the existence of stationary pulses and thus opening way for them to be 
stable, was proposed and studied by means of analytical perturbative methods in Ref. ||, and 
then tested by direct simulations in Ref. [10|. In this system, the CGL equation is linearly 
coupled to an additional dissipative equation, which is a linear one in the most fundamental 
and physically relevant version of the model [10, 11]. In the context of optical fibers, the system 
may be realized as a dual-core fiber, in which an active core carries the linear gain, filtering, 
temporal dispersion, and Kerr nonlinearity, while the parallel-coupled core is lossy, its intrinsic 
nonlinearity, dispersion, and filtering being negligible. It has recently been shown [|l2| that this 
model may describe transmission of fully stable optical solitary pulses with an internal chirp 
(intrinsic phase structure) in an indefinitely long fiber-optic communication link. 

In this paper, we study the stability and collisions of chirped solitary pulses in a system of 
nonlinearly coupled CGL equations stabilized by means of the aforementioned scheme, i.e., by 
linearly coupling each CGL equation to its own linear dissipative counterpart. In its simplest 
versions that are considered in this work, the model includes two or three nonlinearly coupled 
stabilized subsystems, which is a prototype of a stabilized WDM system. The most interesting 
issues, on which we focus in this work, are the stability of SPs in this system and collisions 
between them (including a simultaneous collision between three pulses belonging to three chan- 
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nels). Besides the obvious relevance to optical telecommunications, the obtained results are of 
interest in their own right, demonstrating a new type of stable traveling pulses and collisions 
between them in a generalized (multicomponent) CGL system. 

The paper is structured as follows. In section 2, we give a detailed formulation of the model 
with two channels, and exact solutions for SPs in each channel. In section 3, we develop an 
analytical approach to the collision problem, based on the perturbation theory. In particular, a 
prediction for position shifts of the pulses in the case of a quasi-elastic collision is obtained in a 
fully analytical form. In section 4, results of systematic direct simulations of the collisions are 
displayed for both inelastic and quasi-elastic cases; in the latter case, the analytical predictions 
are found to be in good agreement with the numerical results (in a parametric region where the 
perturbation theory applies). In the same section, a region in the model's parametric space is 
identified, in which the pulses are fully stable, i.e., against both small perturbations and mutual 
collisions. A generalization for a three-channel model is briefly considered in section 5, with a 
conclusion that the pulses are also stable against simultaneous collision of three of them. The 
paper is concluded by section 6. 



2 The two-channel model and exact solutions for the pulses 
2.1 The model 

The simplest version of the model describes the propagation of two waves, u and v, carried by 
two different wavelengths in the active core of a dual-core optical fiber. The waves interact 
through the cross-phase modulation induced by the Kerr effect in the active core. The fields 
u and v are assumed to be linearly coupled to two other fields, <j) and ip respectively, which 
propagate in the passive core placed parallel to the active one. In fact, it is not necessary to 
assume that all the long fiber-optic link has a dual-core structure; instead, it is sufficient to 
have short segments of the dual-core fiber periodically installed into the link. Then, in the same 
uniformly-distributed approximation which was mentioned above in relation to amplification 
and filtering, we may consider an effectively homogeneous dual-core fiber link. In fact, this 
approximation was already applied to the single-channel dual-core model in Ref. [14]. 

Thus, the model is based on the following system of normalized equations governing the 
propagation of the four above-mentioned electromagnetic waves in the two linearly coupled fiber 
cores: 

i(u z + cu t ) + (^D -i S ju tt -iu+ (\u\ 2 + a\v\ 2 )u = K<p, (1) 

% (v z - cv t ) + (\d ~ i) v t t ~ iv + (\v\ 2 + a\u\ 2 ) v = Kip, (2) 

i(& + <y-0t) + »I> = Ku, (3) 
»(^-<J-^) + *T^ = Kv, (4) 

where the complex electric field envelopes u, v, 4>, and tp are functions of the propagation 
distance z and retarded time t which are defined in the usual way M . In the active core, which 
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is equipped with the amplification and filtering, the fields u and v obey Eqs. (||) and @, that 
include the nonlinear self-phase modulation, linear gain, and effective filtering (the latter term 
is formally tantamount to diffusion in the i-space). Coefficients in front of the terms in Eqs. 
dU) and ^ accounting for these three basic effects are all normalized to be equal to 1. On the 
other hand, the dispersion coefficient D is explicitly present in the equations, D > and D < 
corresponding, respectively, to anomalous and normal dispersion. The field envelopes u and v 
in the active core are nonlinearly coupled to each other, interacting through the cross-phase 
modulation (XPM) induced by the Kerr effect, which gives rise to the XPM coupling coefficient 
a = 2 in Eqs. (|) and @. The fields u and v are linearly coupled, through the corresponding 
coefficient K, to their counterparts <j) and tp in the linear dissipative core, which is characterized 
by a loss coefficient V. Finally, the parameters c and 8 in Eqs. (|l])-(£|) stand for the inverse- 
group-velocity differences between the co-propagating waves in the active and passive (lossy) 
cores. 

A numerical value of c (assuming that 5 = c) will play an important role below. To estimate 
the value that is of practical interest, we note that the difference in the inverse group velocity is 
simply related to the frequency separation Aw between the channels, c = (3Au, where (3 = k" 
is the standard physical dispersion coefficient ||, |6| . The frequency difference can be expressed 
in terms of the wavelength separation AA, Alu = — (2-/rnco/A 2 ) AA, where A is the wavelength 
proper, cq is the light velocity in vacuum, and n is the refractive index. For applications, the 
case of interest is the one with the pulse's temporal width ~ 10 ps |2[ |6j, and the wavelength 
separation between 1 nm and 0.1 nm |?J. Using these values, and typical values of the dimensional 
parameters for which the present model provides for the best stability of the pulses (see Eqs. 



(|3q) and (|37|) below), an estimate similar to that presented, e.g., in Ref. |15[ shows that relevant 



dimensionless values of c belong to an interval 

c ~ 20 - 200 . (5) 
Numerical results will be presented for this region of the values of c (see Figs. 6 and 7 below). 
2.2 The linear spectrum 

Before proceeding to the analysis of the full system, it is relevant to consider its linear spectrum 
in the dissipationless limit, i.e., in the case when the gain, filtering, and loss terms are dropped 
in Eqs. (||)-(|^). In the linear limit, the two subsystems (n, (ft) and (v,ip) are decoupled, and, 
looking for the solution to the linearized equations in the ordinary form ~ exp (ikz — icut), one 
arrives at the following dispersion relations between the propagation distance k and frequency 
to: 



1 

Id 



5f-D 2 f±\l {c-5Y-D 2 e +16D 2 K 2 }, (6) 



\2 
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where 

q = k^5-u, £ = u±(c-5)/D, (7) 
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In the definitions (|7|) of the shifted propagation constant and frequency, the upper and lower 
signs pertain, respectively, to the (it, <p) and (v, ip) subsystems, while in the dispersion relation 
@ the two different signs yield two different branches of the dispersion curve, see Fig. 1. 
It follows from Eq. @ that the spectrum shown in Fig. 1 always has a gap, 

< - 4Dq <^(c-5) 4 + 16D 2 K 2 - (c - 5) 2 , (8) 

inside which, following the general principles one may expect the existence of a family of gap 
solitons (note that the full dissipationless version of the present model, including the nonlinear 
terms in Eqs. (|l]) and is definitely nonintegrable, therefore "solitons" are meant here simply 
as solitary waves). However, the objective of this work is not to study that possible family, but 
to focus on the search for stable pulses in the full model, including the gain, filtering and loss, 
which is much more relevant to applications. 



2.3 Solitary-pulse solutions 

If the field is launched into one channel only, the system (|l])-(Q) reduces to a subsystem of 
Eqs. (H) and (||), or @ and @, each including a CGL equation linearly coupled to its linear 
dissipative counterpart. If, additionally, there is no difference in the group- velocity mismatch in 
the active and passive core (i.e., c = 5), then, in the reference frame moving with the common 
group velocity, each subsystem (|l]), (||) or (|2|), (|4|), decoupled from the other one, possesses its 
own pair of exact analytical solutions for chirped SPs, which were actually found in Ref. 11]. 
In the reference frame moving with the common inverse velocity c = 5, the exact solution takes 
the form 



u = Aex.p(ikz) [sech(r/t)] 1+ * M , (9) 
cp = AK{iT - k)~ l e^{ikz) [sech(r/t)] 1+ ^ , (10) 

where the definition of the retarded time t is adjusted to the above-mentioned moving reference 
frame, and the chirp coefficient is 

fi = ~D + ]y/32 + 9D 2 . (11) 
4 4 

The SP's inverse width r] and squared amplitude A 2 (which is the peak power, in the application 
to optics) are expressed in terms of the wavenumber k, 



2 2fc(r-l) 

rf 



k (Dfj, — 2) + T [(1 - n 2 ) D + 4/x] ' 



* (3 + ^ 2 ),„A (121 



and, finally, the wavenumber itself is determined by a cubic equation, 
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(pD - 2) (k 2 - l)k + ((1 - fi 2 ) D + 4/j) (k 2 — T {K — V 2 )) = 0. (13) 

Clearly, physical solutions of Eq. ( |l3| ) are those which yield a real wavenumber k and r/ 2 > 0. 
Physical solutions usually exist in pairs, and only the one with a larger value of the peak power 



may be stable. As it was demonstrated in Refs. ]10| , 11, [l^j , the SP solution with the larger 
amplitude is indeed stable in a fairly vast region in the parameter space (r, K, D). These pulses 
definitely remain stable in the framework of the full system (|]) - @. Indeed, because the two 
above-mentioned subsystems ([[]), (|3|) and (Q), (EJ) are coupled solely by the nonlinear XPM 
terms, the only additional stability condition for a pulse belonging to either subsystem is the 
linear stability of the zero solution in the mate subsystem, which is always the first condition 
imposed on parameters of the eligible model. 

Note that the exact SP solution displayed above can be extended to the full system of the 
four equations (|l|)-(||) if the group- velocity differences vanish, i.e., c = 5 = : the expressions 
©"(0) then yield a solution to the system of the four equations after the transformation 
u, v — > (u,v)/ v3, 4>, if} ~ > (4>i VO / V^, and setting u = v and 4> = i/j. 

In the most general case, c ^ <5, no exact solution for SPs is available, but pulses can be 
found numerically, see below. In any case, pulses generated by the decoupled subsystems ([[]), 
© and ©, (f§ move at different velocities, hence they may collide. The strong XPM-induced 
nonlinear coupling between the channels, together with the dissipative character of the system 
(|l])-(f|), may give rise to complex dynamical behavior as a result of the collisions. 

Our objective in this work is to study in detail collisions between SPs in the system (|l])- 
(Q) and their stability. Note, in particular, that in the case when the inverse-group-velocity 
differences in the active and passive cores are large and nearly equal, i.e., c ~ J > 1, the 
subsystems (|l]), @ and @, (||) nearly decouple, therefore quasi-elastic collisions are expected 
in this case, while at smaller values of the group-velocity mismatch collisions may be strongly 
inelastic. These expectations are corroborated by numerical simulations which are displayed 
below. 



3 An analytical approach to collisions between solitary pulses 
3.1 The perturbation theory 

In the cases of practical interest to fiber-optic telecommunications, the model is far from any 
exactly integrable limit, therefore only direct numerical simulations of collisions between pulses 
(and of their stability), results of which will be summarized in the next section, are really 
relevant. Nevertheless, some qualitative insight into the collision problem can be gained from an 
analytical approach, assuming that pulses may be approximated as quasi-solitons. Within the 
framework of such an approach in its most general possible form, each pulse, in its own reference 
frame (in which an exact solution is given by Eqs. (||) - (13), assuming that 5 = c), is taken as 



A f(r](t - T)) exp(ikz - icot) , (14) 
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= $ g(r](t - T)) exp(ikz - icri) , (15) 
« = Af(rj(t + T))exp(ikz + iut), (16) 
^ = $g(rj(t + T))exp(ikz + iujt). (17) 

Here f(rjt) and g{riT) are (generally speaking, complex) functions accounting for a particular 
shape of the unperturbed pulses, rj being their inverse temporal width, A and <3? are amplitudes of 
their two components, and ±o> and ±T are shifts of the pulses' central frequencies and temporal 
positions due to the interaction between them (we consider the interaction between identical 
pulses, hence the symmetry between the expressions 

(ED, (HI) and (H), (U)). Each component 

of the pulse has its own effective mass, for instance, 

/+oo /»+oo 
\f(x)\ 2 dx, M = $V X / \g(x)\ 2 dx (18) 
-oo J — oo 

(note that, in the absence of losses and gain, M u + M^, as well as M v + M^, are the conserved 
optical energies in each subsystem). For the pulse given by the solution @ - ([l3|), one finds 

M M = 2r 7 - 1 A 2 , = M U K 2 / (k 2 + T 2 ) . (19) 

The XPM-induced coupling between the two subsystems gives rise to a potential force of 
attraction between the pulses, which can be calculated by means of well-known methods (see, 
e.g., Refs. |15|, |H|), provided that XPM may be treated as a small perturbation (conditions for 
applicability of this assumption will be considered below). It is also known that the filtering 
term in Eqs. (|l|) and @ gives rise to an effective friction force, which, in the most general case, 
can be evaluated and combined with the potential force by means of the balance equation for 
the pulse's momentum (as it was done, for instance, in Refs. |I7| ). As a result, one arrives at 
evolution equations for the soliton's position and frequency shifts in the following general form, 

£ = (20) 

duJ 2 M u A 4 , 

where the prime stands for the derivative, and the friction coefficient and interaction potential 



arc 



J^\df(x)/dx\ 2 dx 

K - i^\m\ 2 dx ' (22) 

r+oo 

U(y) = / |/(x-i/)| 2 |/(x + y)| 2 di (23) 



oo 



(the XPM coefficient a was set equal to its physical value 2). The additional term cz in the 
argument of the potential in Eq. (|2l|) is generated by the group-velocity difference between 
the two channels, and the ratio of the masses in the friction term on the right-hand side of 



7 



Eq. (^l]) appears since the friction force acts only on the ^-component of the pulse, but not 
on its ^-component. Note that these general equations are also valid in the case of dispersion 
management, when D is not a constant, but a function of z [17]. 



For the pulses with the shape given by Eqs. (||) - (|i3[) , one can find, from the expressions 
|) and (||), that 



1 + 



k 2 + T 2 



U(y) = 4 ^cosh(2^)-sinh(2,) 
sinh (2y) 



It is worthy to mention that the expression (25) contains no singularity at y — > 0. 



To predict results of the collisions in the general case, the nonlinear non-autonomous (in- 
dependent) ODEs (ordinary differential equations) (|20|) and (^) with the effective potential (25) 



must be solved numerically. In view of the complexity of this ODE system and its approximate 
character, it makes sense to focus, instead, on direct simulations of the underlying PDEs (partial 
differential equations) (|l]) - (||, which will be done below. Nevertheless, some results can be 
obtained directly from ODEs (|20|) and (pl|). In particular, the most essential effect observed 
in direct simulations of the underlying PDEs is an inelastic outcome of the collision (merger or 
complete decay of the pulses), provided that the group- velocity difference c is below a certain 
critical (threshold) value c cr . This value may be, very roughly, estimated as that at which the 
friction and potential forces in Eq. ([y]) are comparable, that yields 

\n\A 2 

lDl 1 (26) 



^(1 + M 2 )' 

To obtain this estimate, it was set that M u j (M u + M^) ~ 1, which is true in the cases considered 



below, the expression (24) for the friction coefficient was used, and it was naturally assumed 
that, for a nontrivial collision, the maximum value of the frequency shift lo is on the order of 
Dc, see Eq. (pc|). Below, it will be seen that the crude estimate (|26|) helps to understand the 
fact that c cr is much smaller for the case of normal dispersion than for pulses propagating under 
anomalous dispersion. 

3.2 Collision-induced position shifts of the pulses 

The ODEs (||) and @) can be used to obtain quantitative results in the limiting case of large 
c, so that 

dT 



iz «c. (27) 

In fact, this is the case when the interaction of the pulses due to XPM may be treated as a 
small perturbation, and all the above approach is strictly valid. In this case, the term rfF in the 
argument of the potential in Eq. ( |2l| ) may be omitted, hence the equation immediately becomes 
linear. Upon Substitution of uj = —D~ 1 dT/dz from Eq. ( p0[ ) into Eq. ([H]) and integrating once, 
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it reduces to the following first-order linear equation: 

^ + K 2(,r) = |# [ /(^) (28) 

a \ cr]Z) c 2c z 

(it is more natural to consider, as final dynamical variables, the renormalized temporal shift r/T 
and propagation distance crjz). To obtain Eq. (|28|), it was again assumed that M u j (M u + MA ~ 
1, which will be the confirmed below, and it was substituted M u = 2 A 2 /rj, as per Eq. (|l9|). 

Equation (^) can be further simplified if, in addition to the condition (|?]), the group- 
velocity difference between the channels is large enough in comparison with an effective friction 
force, so that 

K7] <C c. (29) 



In fact, this condition turns out to be less restrictive than the one (27), see below. Neglecting 
the friction term, Eq. ( |28|) takes the form 

d(riT) A 2 D , . 

^ = -» u m (30) 

Equations (|28| ) and (|30| ) show that the collision is elastic in the present the inverse- 

group-velocity shift dT/dz is zero at z = ±oo, i.e., both before and after the collision. Never- 
theless, the result of the collision is not trivial. Indeed, Eq. (30) can be used to evaluate an 
important characteristic of the elastic collision, viz., a residual temporal (position) shift of the 
pulse, 

A 2 D r + °° 

r]AT = r ] [T(z = +oo)-T(z = -oo)] = — r U(x)dx. (31) 

^ c J -oo 

This shift is important as it gives rise to the collision-induced temporal jitter of the pulses, see, 
e.g., Refs. [15] and [16]. In particular, for the potential ( p5| ) one has U(x)dx = 2, hence 

f] AT = A 2 D/c 2 . (32) 

This analytical prediction will be compared below with results of direct simulations. 

To conclude the analytical consideration, we note that the general condition (|27]) takes a 
simple form in terms of AT. Indeed, in the present case the characteristic collision distance Az 
is determined by the pulse's temporal width 1/rj, so that Az ~ 1/ (rjc), and dT/dz may then be 
estimated as AT/ Az ~ rjcAT. Inserting this into Eq. (|27|), one arrives at a simple result: 



»7|AT|«1, (33) 

which means that the linear equation ( ^8[ ) applies to the description of collisions between pulses 
if the resulting normalized temporal shift of the pulse is small. 

4 Numerical analysis of collisions and stability of pulses 
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4.1 The approach to the problem 



We have employed the split-step Fourier algorithm to solve Eqs. ©-(H) numerically, using, as 
initial conditions, a superposition of separated waveforms @ and (|To|), which yield exact SP 
solutions for the two decoupled subsystems (0), @ and f^), @). Thus, the following initial 
configurations are used: 



u(0,t) = A [sech (77 (t - r))] 1+ ^ , u(0,t) = ^[sech(r ? (t + T))] 1+ ^, (34) 

^(°'*) = ^T^[^ch(r,(t-T))} 1+ ^, ^(0,i) = ^- I [sech(r ? (t + T))] 1+ ^, (35) 

which incorporate an initial temporal separation 2T between the pulses. We will then be able to 
study both incomplete and complete collisions, corresponding to T = and T / respectively 
(i.e., collisions between initially overlapped and separated SPs, see, e.g., Refs. and |[(| for 
the discussion of relative importance of both types of the collisions). 

As for the choice of parameters, in most cases we have used the values T = 5 and K = 4, 
which are located almost in the center of the stability domain of the exact SP solution to the 
decoupled subsystems (|l]), (||) and (§), @ Jl^]. For the dispersion parameter D we have chosen 
the values D = ±18, in the anomalous- and normal-dispersion regimes, respectively. These two 
values of D actually correspond to the carrier wavelength near the zero-dispersion point in a 
dispersion-shifted fiber 0] (recall that Eqs. (|l]) and (||) are normalized so that the effective 



filtering coefficients in them are set equal to 1). As was shown in Ref. [12], for a typical physically 
relevant value of the filtering, the corresponding values \D\ are indeed close to 18. Also, this 
value of D gives rise to the best stability characteristics for SPs in the single channel model. 
Using these values for T, K, and D, the other parameters of the exact-SP solution can be found 
from Eqs. (]9|)-(|l3|): in the case of anomalous dispersion (D = +18), 

\i = 0.074, 7] = 1.57, k = 23.35, A 2 = 44.56, (36) 
and for the normal dispersion (D = —18), 

fi = 27.074, r] = 0.058, k = 23.35, A 2 = 22.37. (37) 

In the anomalous-dispersion regime, the pulses are much narrower, and (quite naturally) have a 
much smaller chirp than their counterparts existing in the case of the normal dispersion. Note 
that, for these values of the parameters, 0.028 • M u according to Eqs. (|l~9|), i.e., the 

mass of the passive-core component of the pulse is negligible in comparison with its active-core 
component's mass. 

Our first objective is to study in detail all possible outcomes of collisions of stable moving 
pulses. In the numerical simulations we have found that, depending on the value of the inverse 
group-velocity mismatch c in the active core, three different outcomes of the collisions occur: 
(a) both SPs perish ("decay"); (b) only one SP survives the collision, while the other one is 
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destroyed (this outcome may also be considered as a merger of two pulses into one); (c) the 
pulses undergo a quasi-elastic collision, so that both reappear unscathed after the collision. 

Formation of a true stable bound state of two solitons as a result of the collision has never 
been observed in the simulations. However, it will be shown below that, in some special cases 
(see Fig. 5(a)), a metastable bound state is observed, which exists over a very long propagation 
distance, but finally collapses into a single pulse. 

As far as the above-mentioned outcomes (a) and (b) are concerned, it is important to mention 
that, in most cases (apart from the exception corresponding to the formation of the metastable 
bound state, which will be specially considered below), we have found that the corresponding 
collision distances are quite short, z co \\ < 60. Thus, in these cases, after passing the short 
collision distance, there remains, at most, one pulse. Obviously, the only outcome acceptable for 
applications to the optical telecommunications is a quasi-elastic collision, when both SPs restore 
their shapes after the interaction. 

Detailed results obtained for the incomplete and complete collisions in both normal- and 
anomalous-dispersion regimes, as well as results for the stability of isolated moving pulses, are 
summarized below. In all the simulations, the XPM coupling coefficient in Eqs. ([!]) and @) was 
set equal to its physical value, a = 2. 



4.2 Incomplete Collisions 

The results for incomplete collisions (T = in Eqs. (34) and (|35|)) in the case c = 5 are 
summarized in Table 1, where the three above-mentioned possible outcomes, namely "Decay", 
"Merger" and "Elastic", are indicated. As is shown, for both the normal- and anomalous- 
dispersion regimes, there exists a critical (i.e., minimum) value c cr of the velocity c, above which 
the collision is always elastic. More importantly, the value of c cr in the case of the normal 
dispersion is much smaller, by a factor ~ 5, than that for the anomalous-dispersion regime. As 
the smaller critical velocity difference between the channels makes it possible to have a denser 
WDM system, this result shows that the normal-dispersion regime may have an advantage over 
the more traditional, from the viewpoint of the soliton transmission [^j, anomalous-dispersion 
regime. On the other hand, an advantage of the latter regime is that, inside a given channel, 
the pulses forming a data-carrying stream may be packed with a higher density, as their width 
is much smaller according to Eqs. (^) and (]37|) . In fact, the best approach to the enhancement 
of the bit-rate of the fiber-optic telecommunication link is to use the channels in both the 
anomalous- and normal-dispersion bands. 

The fact that the critical value c cr is much smaller in the normal-dispersion regime can be 
explained by the crude estimate ( p6|) obtained above on the basis of the analytical consideration. 
Indeed, the ratio of the values which the expression (26) takes for the parameters ( p6| ) and (37) 
corresponding to the anomalous and normal regimes is very small (however, the estimate is too 
crude for a detailed quantitative comparison with the numerical results). 

Typical examples of the three different outcomes of the incomplete collisions are shown 
in Figs. 2(a)-(c) (for D = —18) and Figs. 3(a)-2(c) (for D = +18) in the form of contour 
plots. As is readily observed, SPs in the normal-dispersion regime are indeed much broader 
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than those in the case of the anomalous dispersion, in accordance with what is predicted by 
the analytical solution, see Eqs. (|3^ ) -(|37| ) . In order to test the sensitivity of c cr to a possible 
mismatch of the interchannel inverse-group-velocity differences 5 and c between the two cores, 
we have also simulated incomplete collisions in the case 5 7^ c, viz., 5 = 1.1c and 5 = 2c. The 
results are summarized in Tables 2 and 3, respectively. It is found that, in both normal- and 
anomalous-dispersion regimes, a mismatch in the values of the inverse-group- velocity differences 
in the active and passive cores results in an increase of the critical value c cr , which was again 
found to be much larger in the case of the anomalous dispersion. As was mentioned above, 
for the applications it is necessary to make the critical value c cr as small as possible. The 
results presented here clearly show that the optimum will be attained when 5 = c, i.e., when 
the group- velocity differences between the channels are the same in the active and passive cores, 
which is not difficult to understand in qualitative terms. Indeed, a group-velocity mismatch 
between the two cores makes it necessary for the main component (in the active core) to "drag" 
its counterpart in the passive core, which inevitably generates additional losses through the 
filtering term, thus enhancing inelasticity of collisions between the pulses, cf. the perturbative 
treatment of the collision in the previous section. 

4.3 Complete Collisions 

4.3.1 Inelastic and elastic collisions 

The results of the simulations for complete collisions (i.e., for a case of a sufficiently large 
initial separation between the colliding pulses) in the case 5 = c are summarized in Table 4. 
In the simulations, the initial temporal separation between the solitons is taken to be equal, 
approximately, to five pulse widths, i.e., r[I = 2.34 (see Eqs. (||) and (|To|)). According to Eqs. 
(|36| ) and (^), this choice implies T = 40 and T = 1.5, in the regions of normal and anomalous 
dispersion, respectively. 

The results shown in Table 4 suggest that, contrary to the case of incomplete collisions, in 
the normal-dispersion regime the three outcomes, "decay", "merger" and "elastic", alternate 
with the increase of c, up to the threshold value c cr = 13.92, past which only elastic collisions 
take place. Notice, in particular, a small interval, 7.85 < c < 8.09, where the pulses undergo 
elastic collisions, which is found between the regions where at least one pulse disappears after 
the collisions. Typical examples of the "merger", "decay" and "elastic" outcomes for D = — 18 
are shown in Figs. 4(a) through 4(c). Notice that the collision distance in all the cases does not 
exceed z ~ 5. 

On the other hand, in the anomalous dispersion regime, outcomes of complete collisions 
resemble what was observed in the case of the incomplete collisions: there are only three, 
relatively broad intervals of c, where the outcomes are "merger", "decay" and "elastic", with 
a value of c cr slightly smaller than that in the normal-dispersion regime. However, in the case 
of D = +18 it is important to note that, in the "merger" interval (0 < c < 3.33), the pulses 
demonstrate a behavior which is remarkably different from that observed in the case of the 
"decay" or "merger" outcomes of incomplete or complete collisions in all the other cases. In 
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this interval, the pulses stick together and propagate in such a quasi-bound state over a long 
distance (up to z ~ 3000), after which one of the pulses eventually decays. This behavior is 
demonstrated in Fig. 5(a) for c = 2, where the colliding SPs form the locked configuration 
immediately (see the inset in Fig. 5(a)), and then they propagate, keeping this configuration, 
up to z ~ 3000, where the merger eventually takes place through the destruction of one of the 
pulses. Examples of other outcomes of the collision, namely, decay of both pulses and their 
elastic collision, are shown in Figs. 5(b) and 5(c), respectively, with the collision distance being 
very small, z ~ 0.1. 



4.3.2 Position shifts of the pulses in the case of elastic collisions 

A significant post-collision effect, in the case of the complete elastic collisions between the SPs, 
is a temporal shift AT. The shift is apparent, for instance, in Figs. 4(c) and 5(c). We have 
performed simulations to obtain the normalized temporal shift r/AT as a function of the inverse- 
group-velocity difference c = S, for different initial separations between SPs. The results are 
shown in Fig. 6 (for D = —18) and Fig. 7 (for D = +18); note that these figures display the 
region of the values of c which is relevant to the applications, according to the estimate (||). 
Each curve starts from the maximum value of the temporal shift corresponding to the critical 
value c = c cr . As is seen, the smallest value of the temporal shift is attained in the case of 
the normal dispersion. Also, it is observed that a larger initial separation r]T leads to smaller 
values of the temporal shift in the case of the normal dispersion, while the opposite holds in the 
anomalous-dispersion case. 

It is quite relevant to compare the numerical results displayed in Figs. 6 and 7 with the 
analytical prediction (^), which should be relevant under the conditions (j33|) and (p9|). First 



of all, the substitution of the expression (32) into the condition (|33D leads to an inequality 
c 2 S> A 2 |D|. The values of the parameters being taken as per Eqs. (13(f ) and (|37|), we conclude 
that this inequality holds in the parts of Figs. 6 and 7 where irregular oscillations go over into 
a systematic decay of rjAT with the increase of c. 

Next, using the expression (p3) for the effective friction coefficient which appears in the 



second condition (29), and the same values of the parameters from Eqs. ( |36| ) and ([3?]), one sees 
that this condition readily holds in the case of the anomalous dispersion, shown in Fig. 7, and in 
the case of normal dispersion (Fig. 6) it holds for large values of c, where the above-mentioned 
systematic decay of i/AT takes place. 

Thus, one may expect that the analytical result ( |32[ ) may be correct for sufficiently large 
values of c in both cases. The inspection of Figs. 6 and 7 corroborates this expectation: despite 
a considerable scatter of the values of the normalized shift, depending on the initial value of the 
temporal separation between the pulses, the analytically predicted dependence (|32|) not only 
qualitatively agrees with the numerical results for large c, but also, as one can readily check, 
numerical values of the shift, as predicted analytically and found from the simulations, are fairly 
close. 
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4.4 The full-stability region for the solitary pulses 

The results presented above were restricted to the fixed values of the loss and coupling coeffi- 
cients, r = 5 and K = 4. It is also important to investigate pulse collisions in the full parameter 
space (r, K, c) , and, in particular, to identify parametric domains where solely elastic complete 
collisions occur. These domains actually represent regions of the full stability of SPs in the 
system of Eqs. (||)-(Q), since collisions are natural finite perturbations in this model, against 
which the pulses must be stable, as well as against infinitesimal perturbations (in this extended 
definition of the stability, we take into regard only complete collisions, which are inherent per- 
turbations, and disregard incomplete collisions, that strongly depend upon particular initial 
conditions). 

The thus defined full-stability domains, found by means of systematic simulations, are dis- 
played in Figs. 8 and 9 for the cases of the normal and anomalous dispersion, respectively, as 
gray regions in the (r, K) parametric plane for different values of the group- velocity parameter 
within < c < 10. Notice that the domains have, roughly speaking, a boomerang-like shape, 
resembling the corresponding domain found earlier in the single-channel stabilized CGL model, 
described by the decoupled subsystems ([[]), (H) or @), (|j) fli2| . More importantly, in both cases 
D = ±18, there exists a minimum value of c necessary for the collisions to be elastic, which is 
significantly lower in the case of the anomalous dispersion, namely c m i n = 2 for D = +18, and 
c m in = 6 for D = —18. It should also be noted that the unshaded triangular region shown in 
Fig. 9 corresponds to the case where the two pulses co-propagate undistorted without actual 
interaction over a very long distance (z rj 4000). 

5 The three-channel model 

The model considered in this work can be extended to explicitly include a larger number of WDM 
channels, i.e., a larger number of the CGL equations, each being coupled to its linear dissipative 
counterpart. The simplest generalization contains three channels, which are described by the 
following system of the coupled CGL equations: 



— iu + ( u 


| 2 + 2\v\ 2 + 2\w\ 2 ) u 


= K<f>, 


(38) 


— iw + [\w 


\ 2 + 2\u\ 2 + 2\v\ 2 ) w 


= K X , 


(39) 


— iv + (\v 


| 2 + 2|u| 2 + 2H 2 )u 


= KiP, 


(40) 




i (<f> z + 2c(j)t) + iT<j> 


= Ku, 


(41) 






= Kw, 


(42) 




i (ip z - 2cip t ) + ^rv> 


= Kv, 


(43) 



where we have assumed that each channel has the same group velocity in the active and passive 
cores, while the relative velocity between adjacent channels is 2c, as in the above model (|l|)-(|4|). 
Although detailed study of the extended model is beyond the scope of this work, we give here 
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an example of the existence of fully stable pulses in the system (|38[)-(^3|), which undergo elastic 
complete collisions with each other. As is shown in Fig. 10 (for r = 5, K = 4, c = 20, and 
r/T = 2.34), the three pulses after traveling a distance of z ~ 1 collide and then restore their 
shapes and propagate undistorted. Thus, we may conjecture that the proposed stabilized scheme 
may be generalized to include a larger number of channels in which case stable pulses may still 
experience elastic collisions. 

6 Conclusion 

In this paper we have studied in detail collisions of chirped solitary pulses in nonlinearly coupled 
cubic Ginzburg-Landau (CGL) equations, each being linked to a stabilizing dissipative linear 
equation. Primarily, the two-channel model was considered. The model may be realized as 
wavelength-separated data-transmission channels co-existing in a nonlinear dual-core optical 
fiber, which contains an active core with gain, and a passive core, where the propagation is 
governed by the linear dissipative equations. Each channel has its components in the active 
and passive cores, with a linear coupling between them. Nonlinear interaction between different 
channels is induced by the cross-phase modulation (XPM), which acts in the active core only. 
Thus, the model describes a WDM multi-channel fiber-optic transmission system, provided 
that the dispersion length of the pulses is essentially larger than the amplification and filtering 
spacing, so that the system may be considered in the approximation which assumes a uniformly 
distributed bandwidth-limited gain in the active core (and a continuous passive core, which, in 
reality, may consist of short segments periodically inserted into the long fiber-optic link, together 
with amplifiers and filters). 

If the CGL subsystems are decoupled, they possess stable chirped pulse solutions, which can 
be found in an exact analytical form, provided that the group-velocity parameter is identical in 
the cubic and linear equations. The XPM-induced nonlinear coupling between the subsystems 
gives rise to interactions when the pulses belonging to the different subsystems collide. By 
means of direct simulations, we have studied incomplete and complete collisions in detail. Three 
different possible outcomes of the collision have been found, in the cases when the dispersion 
in the active core is anomalous or normal: destruction of both pulses, destruction of one of 
them, and a quasi-elastic collision. In the latter case, both pulses reappear unscathed after the 
collision (with some positional shifts), provided that the group- velocity difference between the 
channels exceeds a critical value. As a result, regions in the model's parameter space have been 
identified where the pulses are stable against both small perturbations and mutual (complete) 
collisions. An analytical perturbation theory was developed to predict the positional shifts. The 
analytical results agree well with the numerical ones in the cases when applicability conditions 
for the perturbation theory hold. 

For WDM applications, it is important not only to guarantee the quasi-elastic character of 
the collisions between pulses belonging to different channels, but also to have the critical group- 
velocity difference between adjacent channels, necessary for the elasticity of the collisions, as 
small as possible, so that the wavelength separation between the channels may be minimized. 
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To this end, we have found that, in the case of incomplete collisions, the normal-dispersion regime 
provides an essentially smaller critical velocity, whereas in the case of complete collisions, the 
critical velocities are almost equal for both signs of the dispersion. However, the region in the 
parameter space where complete collisions are always elastic is essentially larger in the case of 
anomalous dispersion, and another advantage of the latter case is that the temporal width of 
the pulses is much smaller. On the other hand, the residual effect of elastic collisions, viz., the 
temporal shift of the pulses, which contributes to the soliton jitter in optical communications, is 
weakest in the case of normal dispersion. Actually, the best solution may be to use the channels 
in both normal-dispersion and anomalous-dispersion bands in the fiber. Finally, we have shown 
that the model may be extended to include more than two WDM channels, giving rise, in 
an appropriate region of the corresponding parameter space, to completely stable, three-pulse 
collisions being quasi-elastic. 
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Figure Captions 



Figure 1. A typical form of the dispersion curves q(x) (the shifted propagation constant and 
frequency are defined in Eqs. ([?])) for the linearized system in the dissipationless approximation. 
The curves are shown for D = +18 (the case of the anomalous dispersion) and c = 5. For normal 
dispersion (D = —18), the dispersion curves are mirror images of those shown in this figure. 

Figure 2. Incomplete collisions of two pulses in the normal-dispersion regime {D = —18). 
(a) Decay of both pulses, with c = 8 = 6. (b) Merger of the pulses, with c = 5 = 9. (c) An 
elastic collision, with c = 5 = 14. 

Figure 3. Incomplete collisions of two pulses in the anomalous-dispersion regime (D = +18). 
(a) Decay of both pulses, with c = 5 = 40. (b) Merger of the pulses, with c = 5 = 2. (c) An 
elastic collision , with c = 5 = 48. 

Figure 4. Complete collisions of two pulses in the normal-dispersion regime (D = —18). (a) 
Merger of the pulses, with c = 5 = 10. (b) Decay of both pulses, with c = 5 = 13. (c) An elastic 
collision, with c = 5 = 14. 

Figure 5. Complete collisions of two pulses in the anomalous-dispersion regime (D = +18). 
(a) The merger, with c = 5 = 2. The two pulses get stuck almost immediately (see the inset 
showing the initial stage of the collision in detail), and then they propagate, keeping this shape 
up to z ~ 3000, where the merger (in fact, destruction of one of the pulses) eventually takes 
place, (b) Decay of both pulses, with c = 5 = 5. (c) An elastic collision, with c = 5 = 48. 

Figure 6. Relative temporal shift vs. the inverse-group- velocity mismatch c = 5 in the case 
of normal dispersion {D = —18). 

Figure 7. Relative temporal shift vs. the inverse-group- velocity mismatch c = 5 in the case 
of anomalous dispersion (D = +18). 

Figure 8. The full-stability region (implying the stability of the solitary pulses against both 
arbitrary infinitesimal perturbations, and against collisions with a pulse moving in the other 
channel) in the (K, T) parametric plane, at different fixed values of the inverse-group- velocity 
difference between the channels, in the case of normal dispersion {D = —18). 

Figure 9. The same as in Fig. 8 for the case of anomalous dispersion (D = +18). 

Figure 10. An elastic complete collision between three solitary pulses in the three-channel 
model, in the case of normal dispersion (D = —18), with c = 5 = 20 and rfT = 2.34. 
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Table Captions 

Table 1. Outcomes of incomplete collisions for 5 = c. 

Table 2. Outcomes of incomplete collisions for 5 = 1.1c. 

Table 3. Outcomes of incomplete collisions for 5 = 2c. 

Table 4. Outcomes of complete collisions for 5 = c and r?T = 4.7. 
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Incomplete Collisions (5 = c) 


Normal Dispersion (D = —18) 


Anomalous Dispersion (D = +18) 


Velocity 


Outcome 


Velocity 


Outcome 


c < 8.35 


Decay 


c < 3.52 


Merger 


8.35 < c < 9.53 


Merger 


3.52 < c < 46.98 


Decay 


c > 9.53 


Elastic 


c > 46.98 


Elastic 



Table 1. Nistazakis, Physical Review E 
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Incomplete Collisions (5 = 1.1c) 


Normal Dispersion (D = —18) 


Anomalous Dispersion (D = +18) 


Velocity 


Outcome 


Velocity 


Outcome 


c < 9.88 


Decay 


c < 2.42 


Merger 


9.88 < c < 10.87 


Merger 


2.42 < c < 59.29 


Decay 


10.87 < c < 14.17 


Decay 


c > 59.29 


Elastic 


c > 14.17 


Elastic 





Table 2. Nistazakis, Physical Review E 
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Incomplete Collisions (5 = 2c) 


Normal Dispersion (D = —18) 


Anomalous Dispersion (D = +18) 


Velocity 


Outcome 


Velocity 


Outcome 


c < 14.26 


Decay 


c < 1.73 


Merger 


c > 14.26 


Elastic 


1.73 < c < 122.17 


Decay 




c > 122.17 


Elastic 



Table 3. Nistazakis, Physical Review E 
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Complete Collisions (5 = c and r]T = 4.7) 


Normal Dispersion (D = —18) 


Anomalous Dispersion (D = +18) 


Velocity 


Outcome 


Velocity 


Outcome 


c < 1.88 


Decay 


c < 3.33 


Merger 


1.88 < c < 2.56 


Merger 


3.33 < c < 12.82 


Decay 


2.56 < c < 7.85 


Decay 


c > 12.82 


Elastic 


7.85 < c < 8.09 


Elastic 




8.09 < c < 11.71 


Merger 




11.71 < c < 13.92 


Decay 




c > 13.92 


Elastic 





Table 4. Nistazakis, Physical Review E 
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